On the extended structure of the Phoenix dwarf galaxy 

Sebastian L. Hidalgo^ 
shidalgoOiac . es 

and 

Antonio Aparicio^'^ 
antapajSiac . es 
and 

David Martmez-Delgado^ 
ddelgado@iac . es 

and 

Carme Gallart^ 
carmeSiac . es 

Received ; accepted 



^Instituto de Astroffsica de Canarias. Via Lactea s/n. E38200 - La Laguna, Tenerife, 
Canary Islands, Spain. 

^Departamento de Astroffsica, Universidad de La Laguna/Instituto de Astroffsica de Ca- 
narias. Vfa Lactea s/n. E38200 - La Laguna, Tenerife, Canary Islands, Spain. 



- 2 - 



ABSTRACT 

We present the star formation history (SFH) and its variations with galacto- 
centric distance for the Local Group dwarf galaxy of Phoenix. They have been 
derived from a (F555W, F814W) color-magnitude diagram (CMD) obtained from 
WFPC2@HST data which reaches the oldest main sequence turnoffs. The lAC- 
star and lAC-pop codes and the MinnlAC suite have been used to obtain the 
star formation rate as a function of time and metallicity, ip{t,z). We find that 
Phoenix has had ongoing but gradually decreasing star formation over nearly a 
Hubble time. The highest level of star formation occurred from the formation of 
the galaxy till 10.5 Gyr ago, when 50% of the total star formation had already 
taken place. From that moment, star formation continues at a significant level 
until 6 Gyr ago (an additional 35% of the stars are formed in this time interval), 
and at a very low level till the present time. The chemical enrichment law shows 
a trend of slowly increasing metallicity as a function of time till 8-6 Gyr ago, 
when Z starts to increase steeply to the current value. 

We have paid particular attention to the study of the variations of the SFH as 
a function of radius. Young stars arc found in the inner region of the galaxy only, 
but intermediate-age and old stars can be found at all galacto centric distances. 
The distribution of mass density in alive stars and its evolution with time has 
been studied. This study shows that star formation started at all galactocentric 
distances in Phoenix at an early epoch. If stars form in situ in Phoenix, the 
star formation onset took place all over the galaxy (up to a distance of about 
400 pc from the center), but preferentially out of center regions. After that, our 
results arc compatible with a scenario in which the star formation region envelope 
slowly shrinks as time goes on, possibly as a natural result of pressure support 
reduction as gas supply diminishes. As a consequence, star formation stopped 
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first (about 7 to 8 Gyr ago) in outer regions and the scale-length of the stellar 
mass density distribution decreased with time. Finally, no traces of a true, old 
halo are apparent in Phoenix either in its stellar age distribution or in the stellar 
mass density distribution, at least out to 0.5 kpc (about 2.5 scale-lengths) from 
the center. 

Subject headings: galaxies: dwarf, galaxies: evolution, galaxies: photometry, galaxies: 
stellar content, galaxies: structure 
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Introduction 



During the 50's decade, Baade found an extended structure of red st ars in a series of 



197l|), it 



deep red plates of the IC 1613 dwarf galaxy. Designated as Baade's sheet ( ISandage 
was interpreted as being formed by Population II stars which extended by a factor of 2 with 
respect to Population I stars. During the last years, new exte nded structures have been foun d 



in sev eral dwa rf galaxies or have been confirmed: IC 1613 (IBattinelli. Demers. fc Artigau 



2007 ). WLM f Minniti fc Ziilstra 



1999 ). Antlia (lAparicio et al 



999), Leo A (Vansevicius et al. 



flAparicio fc Tikhonov 



6822 (IKomivama et al. 



996|), NGC 3109 fiMinniti. Ziilstra. fc Alonso 



19971). P hoenix (IMartmez-Delgado. Gallart. fc Aparicio 



20041). DDO 216 flTikhonov 



2003 



2006ah. DDO 190 



2000h DDO 187 (lAparicio. Tikhonov. fc K arachentsev 



Battinelli. Demers. fc Kunkel 



( iHidalgo. Marm-Franch. fc Aparicio 



2000h . NGC 



20061 ). DDO 165, and DDO 181 



20031 ). to name but a few. Other authors have been 



preformed a systematic study of the se extended structures in a wide sample of galaxies wi th 



the aim to derive global properties 


1 

fiDalcanton & Bernstein 


2002; 


Tikhonov 


2005. 


2006b: 


Seth. Dalcanton. & de Jong 


2005a 


b 


). The nature of these structures is still unclear. It 



is, however, of particular relevance to ascertain whether they are or not analogous of the 
halos of spiral galaxies; i.e. formed only by old stars (older than ~ 10 Gyr) and with a hot 
kinematics. Within the context of the standard A-CDM, hierarchical scenario, large galaxy 
stellar halos of such characteristics are predicted to be formed as a consequence of the 
merging of dwarf building blocks. In principle, such a mechanism could also be at work in 
dwarf galaxies si nce, in the same fr amework, they are expected to be scaled-down versions 



of large galaxies (IMoore et al 



200l|). 



From the observational side, 



Minniti fc Zijlstral (119961 ) and 



Minniti et al. 



(119991 ) 



concluded that the extended structure they found in WLM is composed only by old stars 
and that it is a true halo analogous to those of spiral galaxies. However, from their analysis. 
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based in color-magnitude diagrams (CMDs) that are not deep enough to reach the oldest 
turn-offs, it is not possible to exclude the presence of intermediate-age stars in the external 



regions of W LM. The same is true for other dwarf galaxies studied since then ( lAparicio et al. 



1997 



2000l ). In addition, there are several independent observations that hint to the 



presence of intermediate-age populations in these low surface brightness envelopes of dwarf 
galaxies: i) Carbon stars, which are bright members of the intermediate-age population, 
have been identified in dir r galaxies, out to where young stars have already disappeared 



( iBattinelli fc Demers 



2000 



Demers. Battinelli. fc Artigau 



2006 



Battinelli fc Demers 



20061); 



i i) PNe have been a l so detected at lar ge distances from the center of many local dirrs 



( iMagrini et al. 



2003 



Leisv et al. 



brightness fluctuations. 



Hidalgo et al. 



20051 ) : iii) from population synthesis based on surface 



(120031 ) found indications that the extended regions of 
dwarf galaxies could contain a mixture of intermediate-age and old stars, the average age of 
such mixture possibly increasing with galactocentric distance; iv) in irregular galaxies for 



which CMD reaching the oldest main-sequence turnoff exist, such as the Mage^ 



an intermediate-age populat ion has been found in the outskirts of the galaxy (IGallart et al. 



2004bl : 



Noel fc Gallart 



20071). 



lanic Clouds, 



A stellar 'halo' composed by both old and intermediate-age stars is difficult to reconcile 
with a true halo formed through the standard merging of smaller systems in a A-CDM 
framework: since most merging happens in the early history of the Universe, it is unlikely 
that merging produce significant intermediate-age stellar populations in halos. Alternatively, 
mixed-age 'halos' could be produced if stars form in the inner regions of the galaxy and 
migrate to the external ones as a consequence of velocity d ispersion or heating p roduced by 



Hidalgo et al. 



(120031 ) proposed a 



internal or external interactions, among other mechanisms, 
'shrinking' scenario in which the presence of intermediate-age populations in the outer parts 
of dIrr galaxies could be explained if a progressive shrinking of the star formation region, 
initially extending out to the galaxy outskirts, occurred in the galaxy, possibly produced 
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by the gas progressively thinning out in the outermost regions. IStinson et al.l (120091 ) have 
shown, through SPH simulations, that a contracting star formation pattern results naturally 
as a consequence of gas pressure support reduction as gas supply is consumed in the center 
of the galaxy. Stellar migration from the center to the outer parts also plays a role, even 
though it is minor in general. 

Extended structures require a detailed study on an observational basis in order to 
ascertain their nature. First the stellar age and metallicity distributions as functions of 
the galactocentric distance must be determined. These distributions can be found by using 
deep CMDs (reaching the oldest turn-offs) to derive the star formation history (SFH). For 
galaxies beyond the Milky Way halo, the HST is required for this. Second, the kinematics 
in the inner and outer regions of the galaxies must be determined and compared with 
dynamical model predictions. Ground based, 10-meter class telescopes can reach the 
required accuracy for very nearby galaxies. 

This paper is devoted to the first topic, namely the age and metallicity distributions 
of the stellar population or, in other words, the star formation history ( SFH), as a functi on 



of gal actocentric distance in the Phoenix dwarf galaxy. Discovered by 



Schuster fc West 



19761). it is a low surface brigh tness transition galaxy (dSph/dIrr) with /iB = 24.9 mag/( 



il\2 



( Martfnez-Delgado et al 



( iHeld. Saviane. &: Momany 



1999f). At a distance of ~ 450 kpc from t he Milky Way 



1999 



Holtzman. Smith. &: Grillmair 



20001 ) and ~ 600 kpc form 



M31, Phoenix is close enough for observations to reach the old main sequence turn-offs 
but still far away enough to be unaffected by strong gravitational interaction with the 
Milky Way or M31. Phoenix contains a young stellar populati on with strong evidence o 



a population gradient as seen from gro u nd-based observations ( 



Held 



et al. 



199S 



Ortolani fc Gratton 



Holtzman et al.ll2000i ). IMartinez-Delgado et al 



(119991) and 



1988 



Held et al 



( 119991 ) found that the younger and bluer stars of Phoenix are concentrated in the center of 
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the galaxy. The same authors also found that RGB stars are prese nt throughout the whole 



body of Phoenix. Associated with this stellar population gradient, 



Martmez-Delgado et al 



( 119991 ) found a sharp rotation in the elliptical isopleth of ~ 90° at radius r ~ 115" from the 
center, suggesting the existence of two components: an inner one elongated in the east-west 
direction, which contains all the young stars, and an outer component oriented north-south, 
which lacks young stars. Limited by the depth of their photometry, they were unable to 
determine the nature of the stellar populations in the outer component. Phoenix has no H I 
emission coincident with its main body which suggest that all the gas of the galaxy has been 
removed or exhausted by the star formation processes. There are some H I c l ouds detected 



betw e en 4.5' and 9^ from the galaxy center (jCarignan. Demers. fc Cote 



1997 



2001 



Young et al 



2007 



Irwin fc Tolstoy 



Held 



2002). 



1991 



2002) that seem to be related to Phoenix ( 



Young et al. 



Young fc 



Gallart et al 



(120071 ) have used high sensitivity and resolution 



imaging with the VLA to characterize the morphology of the gas cloud and its mass. The 
cloud, with a mass Mhi — IO^Mq, has a crescent-like shape with the center of curvature in 
the direction of the galaxy, which implies that the gas was likely expelled from the galaxy 
by supernovae winds. These characteristics make Phoenix a case of special interest both 
to characterize in detail its full SFH and its variation with galactocentric radius. This can 
provide important insight on the nature of the extended structures in dwarf galaxies. Main 
properties of Phoenix are summarized in Table [H 

The organization of this paper is as follows: In ^the data, the photometric procedures, 
and the calibration method are presented. In ^the CMD is presented. In ^the distance 
to the galaxy is obtained using the tip of the Red Giant Branch (TRGB). In ^ three former 
candidates to globular clus ters are discussed. Sectio n [6] describes t he completeness test s 



perfo rmed. The lAC-star ( lAparicio &: Gallart 



20091 ) algorithms, together with the suite MinnlAC (see also 



Aparicio fc Hidalgo 



2004J ) and I AC-dod (| Aparicio fc Hidalgo 

2009^ have 



been used to obtain the SFH. It is presented in ^ for the whole covered field and in ^ 
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as a function of galactocentric radius. The evolution of the galaxy structure with time is 
discussed in §3 Finally, the main conclusions derived in this paper are summarized in §1U[ 



Reduction and calibration of data 



The images of Phoenix were obtained in January 2001 using the WFPC2 aboard the 
HST (cycle 9, proposal ID 8769, P.I. A. Aparicio). The WFPC2 has four 800 x 800 pixel 
CCDs; three of them (the Wide Field Camera, WFC) with a pixel size of 0.1" pixel"^ and 
the remaining one (the Planetary Camera, PC) with a pixel size of 0.046" pixel"^. The 
resulting field of view are 150 x 150 (")^ for the WFC and 35 x 35 (")^ for the PC. Images 
were taken in two fields observed in the HST filters F555W and F814W: an inner field 
centered on the main body of the galaxy (observed on 2001 Jan 10 and 11), and another 
located ~ 2.7' from the former (observed on 2001 Jan 16 and 17). 



Images obtained by lSmith. Grillmair. fc Holtzmaru (119971 ) (cycle 6, proposal ID 6798, 
P.I. G. Smith), which mat ch our inner f ield d escribed above, were also used. Total 



integration times including 



Smith et al. 



(119971 ) data, were 12100 s and 14500 s for the 
F555W and F814W filters respectively for the inner field; and 11900 s and 14300 s for the 
F555W and F814W filters respectively for the outer field. The observation log is given in 
Table [21 Figure [1] shows the observed fields overplotted on an image of the galaxy obtained 
with the 100-inch Du Pont telescope at Las Campanas Observatory (Chile). 



El 



E2 



EiNOTE TO EDITOR: PLEASE, PUT TABLE |2] HERE 
^^NOTE TO EDITOR: PLEASE, PUT FIGURE □ HERE 



- 9 - 



Processed images of Phoenix were retrieved from the HST archive. There is a small 
difference in the pixel size of the WFPC2 CCDs. Therefore, flux correction per pixel 
was necessary using an image mask kindly provided by P. B. Stetson. Bad pixels were 
corrected with anoth er image mas k . The process to obtain the photometry follows the 



method described by 



Turner et 



al. 



(1l998l ) and is based on the application of DAOPHOT 



and ALLFRAME (iStetsoniil994j ). Moffatian PSFs of parameter f3 = 1.5 and 20 pixel radius, 



provided by P. B. Stetson, were used. 

Photometry was corrected for the charge transfer efficiency (CTE) effect. The 
correction depends on the position of the stars in the CCD, its temperature, and the 
observation date. CTE correction must be applied to the total m agnitude of the st ars, i.e.. 



Hill et al 



Hill et al. 



(1998) were 



(11998h and the 



after aperture corrections. The aperture corrections proposed by 
used. The difference between the focus of the HST for the data of 
data used here is 4 fim, which corresponds to a difference in the flux of the stars of less 
than 0.5% . This small difference lends support for us to assume the correction apertures of 



Hill et al. 



(119981 ) to estimate the CTE corrections for our photometry. 



The final photometric list was cleaned of non-stellar objects and badly measured 
stars by limiting the a and SHARP values provided by ALLFRAME (a < 0.2, 
-0.5 < SHARP < 0.5). No limits in CHI (also provided by ALLFRAME) value were 
imposed because of the restrictions in a eliminated those stars with large CHI. Figure [2] 
shows the errors provided by ALLFRAME as a function of magnitude. 



E3 



The photometric calibration was carried out using V and I ground based observations 
of Phoenix from the 100-inch Du Pont telescope at Las Campanas Observatory (Chile) 



^^NOTE TO EDITOR: PLEASE, PUT FIGURE [2] HERE 
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( IGallart et al.ll2004bl ). A selection of ~ 1600 bright and isolated stars in Phoenix observed 
from the ground were used as local standard stars. The photometric transformation 
equations were obtained comparing the magnitudes of local standards in V and / 
Johnson-Cousins filters with the magnitudes of the same stars observed in the HST F555W 
and F814W filters. For stars of spectral types 05V to M6V, these f ilters are very close 



Bagget et al. 



2002). 



approximations to the Johnson- Cousins filter V and Ic respectively ( 
The WFPC2 calibration was done independently for each CCD. The final photometric 
calibration equations can be expressed as: 



V-v = Zv + CviV -I) + KvY 
I -I = Zr + Ci-{V -I) + Kj-Y, 

(1) 

where V and / stand for calibrated magnitudes and v and i are the instrumental 
magnitudes. The coefficients Zyj and Cyj are, respectively, the zero point of the 
photometry and the color term. The Kyj coefficient is a term that accounts for the 
dependence of the magnitude with the Y coordinate of the star on the CCD. The correction 
introduced by this term makes the magnitude of a star dimmer with increasing Y (CTE has 
the opposite effect). This could be either a new term in the photometric calibration of the 
WFPC2 or a correction to the CTE. Table [3] lists the values of the coefficients for each CCD 
and filter (column 1). The photometric zero point (Z) is given in column 2. Column 3 gives 
the color term (C). The term K for the dependency on the Y coordinate is given in column 
4. Finally, the photometric transformation error is given in column 5. There are two sources 
contributing to the latter: i) the error of the ground based photometric transformations, 
which is the dominant, and ii) the error of the transformation from the HST to ground based 
photometry. Both have been estimated as the rms of the corresponding data distributions 
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divided by \/N — 2, where is the total number of measured stars in each case. The error 
given in column 5 of Table [3] is the quadratic sum of both. 

E4 



3. The CMD of Phoenix 

The calibrated CMD of Phoenix is shown in Fig. [3l Absolute magnitudes are given 
on the right by using a distance modulus of (m — M)o = 23.09 mag and extinctions 
Ay = 0.062 mag and Aj = 0.036 (see §1]). The CMD shows a clear main sequence (MS) 
and RGB. The brightest stars in the MS are younger than 1 Gyr, while the RGB is formed 
of stars older than 1 — 2 Gyr. Parallel to the blue side of the RGB, some AGB stars can 
be observed. A few stars above the TRGB could be intermediate-age AGB stars. An 
extended Horizontal Branch (HB) is also present. This is an indication of very old, metal 
poor stars. The gap produced in the HB by the RR-Lyrae variables is also visible. Near the 
intersection of the HB with the RGB, the Red Clump (RC) is observed. Finally, some blue 
core helium-burning (Blue-HeB) and red core helium-burning (Red-HeB) stars are seen 
extending to brighter magnitudes above the RC. These stars are younger than ~ 0.5 Gyr. 
In short, the CMD of Phoenix shows stars of all ages. 

E5 



NOTE TO EDITOR: PLEASE, PUT TABLE |3] HERE 
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Distance to Phoenix 



The distance to Phoenix has been estimated 
the prescriptions method by 



Bellazzini et al 



a Sobel filter of kernel [1,2,0,-2,-1] (IMadore &: Freedman 



rom the TRGB magnitude, following 



(l2004f). For our d etermination of the TRGB, 



19951 ) was convolved with the 



luminosity function of the RGB stars redder than {V — I) = 0.8 with magnitudes 
18.5 < / < 21. The luminosity function was obtained using a binning of 0.1 mag with a 
step of 0.05 mag. Figure H] shows the result of the convolution. The TRGB is found to be 
at Itrhr = 19.14 =b 0.02 mag. Us ing extinctions Ay = 0.062 mag and Aj = 0.036 mag 



( ICardelli. Clayton. &: Mathis 



19891), the TRGB corrected of reddening is found to be at 



iTRGBfi 



found by 



19.10 ± 0.02 rnag wh ich is, within the errors, located at the same magnitude 



Holtzman et al. 



(120001 ) . The error was estimated as the half maximum of height 



of the convolution peak. 

The absolu t e mag nitude of the TRGB has been calculated using the results by 



Bellazzini et al. 



(120041 ) and assuming [M/H] = [Fe/H]. The metaUicity [Fe/H] has been 



estimate d using the mean co 



given by lSaviane et al. 



or of the RGB stars at Iq = —3.5 and the empirical calibration 
(l200o[ ). The RGB mean color is {V - /) -3.5,0 = 1.32 ± 0.01, which 
gives [Fe/H] = -1.87 ± 0.06 dex (Z = (2.6 ± 0.3) x 10"^). Using [Fe/H] = -1.87±, 
we obtain an absolute magnitude of the TRGB MJq^^ = —3.99 mag which results in a 
distance modulus of (m — M)o = 23.09 ± 0.10, or a distance d = 415 ± 19 kpc. The errors 
include the uncertainties in the calibration method and the determination of the magnitude 
of the TRGB. 



E6 



The distance obtained here is in full agreement with the distances obtained by 



^^NOTE TO EDITOR: PLEASE, PUT FIGURES S HERE 
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van de Rvdt. Demers. fc Kunkell fll99lh ((m-M ) ^ = 23 .1±0.1), by 



1999h ((m- M 



Holtzman et al. 



Menzies et al. 



23.0 ± 0.1), by 



Held et al. 



Martmez-Delgado et al. 



fll999h ((m - M)o = 23.21 ± 0.08), by 



(120001) ((m — M)o = 23.11), and with the most recent distance estimation by 
( I2OO8I ) using infrared photometry and variable stars ((m — M)o = 23.1 ± 0.2). 
The obtained metallicit y ([Fe/H] = —1.8 7 ±0.06 dex) is also in good agreement 



with those obtained 



Ho 



by 



tzman et al. 



by 



Held et al. 



fll999h ([Fe/H] 



-1.81 ± 0.10 dex), and by 



(120001) f[Fe/H] = —1.81 dex), but it is lower than the metallicity obtained 



Martmez-Delgado et al. 



(I1999I ) ([Fe/H] = —1.37 dex) using ground based observations. 



This discrepancy for the metallicity coul d be explained by a mistake in the sign of one of 



the color terms in the equations used by 
instrumental magnitudes into the Johnson-Cousins system (see 
affects the color of the RGB. 



Martmez-Delgado et al.l (119991) to transfor m the 



Gallart et al. 



2004a) which 



5. The globular clusters candidates in Phoenix 



Globular clusters (GC) detected in the Local Group s ^ 



low a luminosity distribution in 



I991I ). They are usually selected 



the range of My > -7.1 and 0.7 < {V - I)o < 1.9 (iHarrisI 
from images in which they are unresolved using these criteria as well as the PSF profiles. 
However, these characteristics may also correspond to background galaxies. 



T here are currently six globular clusters candidates in Pho enix identified in previous 



works (ICanterna fc Flower 



1977 



Martmez-Delgado et al. 



19991). We have ide ntified in 



our HST images the three objects proposed by iMartmez-Delgado et al. 



(11999h . They are 



shown in Figure [5l Visu al inspection indicates tha t the objects are background galaxies 



Candidates identified by 



Canterna fc Flower 



(I1977I ) are out of our observed field. 
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6. Completeness tests 



Our analysis o f the p hotometry completeness follows the procedure proposed by 



Aparicio fc Gallart 



(Il995l ). In short, a list of artificial stars is made. The absolute 



magnitudes of these stars are transformed to instrumental magnitudes with the equations 
used to calibrate the instrumental photometry (CTE, aperture corrections, zero points and 
color terms) and assuming the distance modulus and reddening given in §4l The indiyidual 
PSF for each frame is used to inject the artificial stars into the single frames. 

A synthetic CMD wi th 500000 stars was used as catalog of artifi cial stars. It was 



gener ated with lAC-star (lAparicio fc Gallart 



20041 ) using the B aSTI ( 



Pietrinferni et al 



CasteUi fc Kuruczl (12004] ). 



20041 ) stellar evolution library and the bolometric corrections by 
A constant star formation rate (SFR) between and 14 Gyr was assumed. An uniform 
range of metallicity 0.0001 < Z < 0.002 for all ages was used. The criterion in the selection 
of the lAC-star input parameters is that colors, magnitudes, ages, and metallicities of the 
artificial stars cover the ranges shown or reasonably expected for the observed stars. This 
is important in order of properly simulate observational effects on other synthetic CMDs. 
The artificial stars were distributed over the 8 CCD chips, the number of injected stars in 
each field being proportional to the number of observed stars. They were distributed in a 
triangular grid with a separation between star centers of 2 x -Rps/ + 1 pixels, where Rpsf 
stands for the PSF radius. This separation prevents the superposition of artificial stars to 
each other. 



Once the artificial stars were injected, the new synthetic frames were processed 



^mOTE TO EDITOR: PLEASE, PUT FIGURE HERE 
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photometrically in the same way and with the same parameters as the original images. 
The same photometric calibration corrections were applied to the magnitudes of the fake 
recovered stars, with the exception of the CTE. Since this is a process that affects the read 
out of the CCD, it is not expected to affect the artificial stars, so we did not include it in 
the photometric corrections for these stars. 

Figure [6] shows the CMD of Phoenix overplotted on the recovered artificial star CMD. 
The completeness factors are shown on the right Y -cLXlS clS db function of magnitude. It can 
be observed that the recovered artificial star CMD fully covers the observed CMD, fulfilling 
our previous requirement. 

E8 



7. The Star Formation History of Phoenix 



We have obtained th e SFH of Phoenix following the method explained by 



Aparicio fc Hidalgol ( 



lAC-pop (see 



20091). It ma 



ces us e of the suite of routines MinnlAC and the algorithm s 



Aparicio fc Hidalgo 



2009 



for both) and lAC-star (lAparicio fc Gallart 



20M). 



he m ost important characteristics of our method (see below and lAparicio fc Hidalgo 



(120091 ) for a more cor nprehensive d i scussi on) compared with others SFH recover ing 



methods as MATCH (IDolphin 



1997 . 



20021) or Star-FISH (IHarris fc Zaritskvll200lh are: i 



interpolation in both age and metallicity directly on stellar evolution tracks to produce 



synthetic CMD with smooth sampling of 
sets of stellar evolut ion models (currently 



3oth coordina t es; ii) possibility to use s evera^ 



Bertelli et al. 



Girardi et al. 



(11994f ) 



Pietrinferni et al. 



(120041 ) 



(I2OOOI )): iii) the powerful n ainimization techniq ue implemented in lAC-pop 



(which makes use of a genetic algorithm (iCharbonneau 



19951 )). allows us to solve for a huge 



^^NOTE TO EDITOR: PLEASE, PUT FIGURE HERE 
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number of parameters (e.g. well sampled age-metallicity relation); iv) different areas of 
the synthetic and observed CMDs are sampled with different resolution, according to the 
reliability of our theoretical knowledge of these areas; v) each SFH solution is given as the 
average of a number of solutions obtained from multiple samplings of both a) the synthetic 
CMD in terms of simple stellar populations and b) the CMDs in terms of boxes in the 
color-magnitude plane. 

We define the SFH as a distribution function ?/'(t, z)dtdz that can be identified with 
the usual SFR but with dependence both on time and metallicity. An arbitrary tplt, z) 
can be written as a linear combination of the SFRs, ipi, of simple stellar populations, each 
containing stars with ages and metallicities within small intervals: 



^{t,z)= AY.ailpi (2) 

i 

where the SFR is constant within each ■j/'j and A is a scaling constant. 

lAC-pop derives z) of the target object by comparing the star distributions in 
the observational CMD (oCMD) and in a sy nthetic CMD. A detailed description of the 



method is given in iAparicio fc Hidalgo! (l2009l ). In short, the procedure is as follows. A large 



synthetic stellar population with ages and metallicities spanning ranges wider than those 
expected for the target object is computed together with its CMD. We denote this starting 
global CMD as sCMD. This population is divided in a set of simple stellar populations, ipi. 
Each one has an associated CMD, which distribution of stars as a function of color and 
magn itude is parametrized using a combination of grids (see below and 



Aparicio fc Hidalgo 



20091 ). The CMD parametrization of simple population i is denoted as M/, where j refers to 
the j-th grid element. Using this information, the parametrization of the star distribution 
in the synthetic CMD of any SFH that can be written according to Eq. (2), can be obtained 
as = A'Y^^a.iMl . It should be noted that no need of explicit computation of the new 
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synthetic CMD is required for this and that it is enough to work with the Ml counts. A 
similar parametrization, denoted O*, is done for the observational CMD. At this point, the 
Ui coefficients of the SFH of the targ et object are o btained by minimizing a merit function 
comparing M/ and OK Mighell's ( Mighelj 1999 ) is used for this purpose. 

Concerning the software, lAC-star is used for the computation of synthetic populations 
and their CMDs; lAC-pop computes the distributions and search for the solutions, and 
MinnlAC takes care of the several CMD parametrization as well as the solution handling 
and error computations (see below). 

Specifically, the steps followed to obtain the SFH of Phoenix were as follow: 



The BaSTI (jPietrinferni et al.ll2004l ) stellar evolution library was used to compute 
the sCMD. It contains 5 x 10^ stars with a constant SFR between and 13.5 Gyr 
and an uniform range of metallicity 0.0001 < Z < 0.002 for all ages. These limits 
were selected to cover the full expected ranges of stellar ages and metallicities. Six 
fractions of binary stars were considered: 0, 20, 40, 60, 80 and 100 % with mass ratio 
randomly distributed between 0.5 and 1. Bolometric corrections by 



Castelli &: Kurucz 



(120041) were used in all th e case s and the Initial Mass Function (IMF) was that by 



Kroupa. Tout, fc Gilmord (119931 ). 



Observational effects were simulated in all sCMD using the results of the experiments 
described in ^ The simulation is done as follows: for each synthetic star in sCMD 
with magnitude m^, a list of artificial stars with |mj — ms\ < e is selected for 
each filter, being e an input interval, e = 0.2, which is the largest error in the 
photometry, has been used. A single artificial star is randomly selected. If m'^ and 
m'^ are the injected and recovered magnitudes of the selected artificial star in a 
given filter, then ml = rUs + m'- — m'^ will be the magnitude of the synthetic star 
with observational effects simulated. The same is done for both filters. With this 
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procedure, the observational effects are simulated in sCMD star by star. The quality 
of the observational effects simulations is determined by the number of injected stars 
in the completeness tests in comparison with the number of synthetic stars in sCMD. 

Simple stellar populations were selected from sCMD by using a 0.5 Gyr age binning for 
stars with ages — 1 Gyr; 1 Gyr for stars within 1 — 6 Gyr years old, and 2 Gyr for stars 
older than 6 Gyr. The metallicity bin size increase with metallicity: Az = 0.0001 for 
stars with z < 0.0005; Az = 0.0002 for stars with 0.0005 < z < 0.0007; Az = 0.0003 
for stars with 0.0007 <z < 0.001; and Az = 0.0005 for stars with z > 0.001. 



According to lAparicio fc Hidalgo! (120091 ) , an optimum parametrization of the CMD 
stellar distribution for the purpose of deriving the SFH is obtained defining several 
significant regions (or bundles) in the CMD and sampling them with grids which bin 
size is fixed for each bundle. In our case, two bundles have been defined (see Fig. [7]); 
one of them on the MS and sub-giant region and the other one on the reddest part of 
the RGB. The former includes ~ 65% of the observed stars and was sampled by a 
grid with boxes of size 0.015 x 0.4 in color and magnitude. The latter contains only 
~ 100 observed stars and was sampled with boxes of size C x 0.4, with C covering the 
full color interval of the region for each magnitude. 



E9 



To minimize the dependency of the SFH on the selection of the simple stellar 
populations and position of the sampling grids, MinnlAC generates several sets of 
grids to parametrize the CMDs and several sets of simple stellar populations. This 
is done by introducing offsets in the positions of the CMDs sampling boxes and 
by moving the limits b of the age and metallicity bins defining the simple stellar 



^^NOTE TO EDITOR: PLEASE, PUT FIGURE [7] HERE 
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populations to 6 ± 0.25 56 where 6b is the bin size. 

Uncertainties in distance, reddening and photometric zero-point may have significant 
effects on the SFH solution. To overcome this, several offsets have been applied to 
the observational CMD and the SFH has been found for each of them. Offsets were 
applied according to a grid of 5 x 5 points separated by intervals of 0.02 in (V — I) 
and 0.05 in I. 

In summary, for each binary fraction and CMD zero-point offset, the different sets of 
simple stellar populations and boxes defined above give a total of 16 parametrisation 
of the CMDs (8 sets of simple stellar populations x 2 grid offsets). The average of 
the 16 solutions is assumed as the final solution for the corresponding binary fraction 
and CMD zero-set offset, and the rms dispersion is assumed to be the associated error 
(see lAparicio fc Hidalgo 



20091 ). The average, xL' '^^ 16 x is also computed and 
considered an indicator of the goodness of the average solution. The solution for the 
binary fraction and CMD zero-point off-set producing the smallest xL is assumed to 
be the solution for the SFH of Phoenix. 

The final solution has xiv — 1-9 ^"^^ corresponds to 40% of binaries and CMD offsets 
of A(y — I) = 0.03 and A/ = 0.075. It is interesting to note that the solution (and its xD 
does not change significantly for binary fractions in the range of 40% to 80%. Also, the 
/S.{y — I) and A/ found should not be used to correct distance or reddening, since they are 
also accounting for possible errors in the stellar evolution library and in the photometry 
zero-point. 



The population box (IHodge 



19891 ) for the SFH is shown in Fig. [HI The height of the 
curved surface over the age-metallicity plane gives il){t, z). The total mass converted into 
stars per unit area (M© pc~^) in the age-metallicity interval ([tj, U + 5t];[-2j, Zj + 5z]) is the 
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volume below the solution and within that interval: rriij = J^' J^^ ^ ilj{t, z)dtdz. The 
SFH as a function of time ip{t) and metallicity ip{z) only are also given as projections of 

ElO 

For more clarity, Fig. [9] shows three different aspects of the SFH of Phoenix. In the 
upper panel the SFH as a function of time only ^(t) is shown. It is the same represented 
in the right-hand side, vertical panel of Fig. [SI Thin lines correspond to the rms dispersion 
computed as explained above. The middle panel shows the chemical enrichment law (CEL) 
in a logarithmic scale. It is obtained from the averaged metallicity of the solution for each 
age. More precisely, this averaged metallicity is: 



Z{t) 



J^^--^it,z')dz' 



(3) 



where Z^ax = 0.002 is the maximum metallicity considered. The variable represented in 
Fig. [n]is the logarithm of Z{t) in solar units, assuming Zq = 0.0192. Thin lines show the 
interval containing 68% of the former integral. The mean metallicity for stars older than 
~ 2 Gyr obtained from Z{t) is Z = 0.0004 ± 0.0001, which is compatible with the mean 
metallicity of the RGB stars obtained in ^ using other method. Finally, bottom panel 
of Fig. [9] shows the normalized cumulative mass converted into stars obtained by simple 
integration of 4'{t). 

We have obtained the SFH of Phoenix using models from the BaSTI stellar evolution 



library. However, a different set of ste . 



different solution (lAparicio fc Hidalgo 



lar evolution mode 



2009 : 



Gallart et al 



s cou ld produce a significantly 



20051 ) . To check the stability 



of our solution, we have calculated once more the SFH of Phoenix but using the stellar 



'^i^NOTE TO EDITOR: PLEASE, PUT FIGURES [HI [9] AND [10] HERE 
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evolution models from Padua flBertelli et al.lll994l ) as input for lAC-star to compute the 
sCMD. The rest of input parameters are the same as those corresponding to the adopted 
solution, in particular 40% of binaries and observational CMD offsets of A{V — I) = 0.03 
and AJ = 0.075. The solution based on the Padua library shows a faster SFR during the 
initial 2-3 Gyrs and is very similar from that epoch on. Differences are relatively small and 
do not affect the conclusions of this paper. Table H] gives a summary of the mean values 
derived from il){t,z) for the BaSTI and Padua libraries. Column 1 identifies the stellar 
evolution library used. Column 2 gives the integral of il){t, z) for the complete range of age 
and metallicity, that is, the total mass ever converted into stars. Columns 3 to 5 give the 
mean value for ip{t), metallicity, and age. 



Ell 



Our results (Figs. [8] and [9]) indicate that the Phoenix dwarf galaxy has had ongoing 
but gradually decreasing star formation over nearly a Hubble time. The highest level of star 
formation occurred from the formation of the galaxy and till 10.5 Gyr ago, when 50% of the 
total star formation had already taken place. From that moment, star formation continues 
at a significant level until 6 Gyr ago (and additional 35% of the stars are formed in this 



time interval), and at a very low level till the present time (t 



are known to b e ^ 100 Myr old ; 



Holtzman et al 



Martmez-Delgado et a 



mm and 



Young et al 



le youngest Phoenix stars 



1999). These results are in good 



mm who used t he WFPC2 data 



agreement with 

for th e inner field of the galaxy (somewhat shallower dataset in the case of 



Holtzman et al. 



20001 ). Their intermediate-age star formation rate, however, seems to be somewhat higher, 
and in fact, 80% of the stellar mass had formed 3 and X Gyr ago, respectively. The fact 
that their SFH is restricted to the inner field is however consistent with this trend, since we 
also find more intermediate-age star formation toward the center of the galaxy (See 



'^iiNOTE TO EDITOR: PLEASE, PUT TABLES] HERE 
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The CEL shown in Figure [9|, central panel, shows a trend of slowly increasing metallicity 
as a function of time till 8-6 Gyr ago, when Z increases steeply to the current value. From 
this CEL, we will try to extract some information about the physical conditions of the 
chemical enrichment in Phoenix. In Figure [TOl the Phoenix CEL is compared with several 
simple metallicity law models. It is apparent that none of the models adequately reproduce 
the Phoenix CEL. However, a few constraints may still be set. Infall models are the most 
clearly departing from the Phoenix CEL along the full history of the galaxy. We can 
therefore speculate that such process has likely not been relevant in the Phoenix history. 
A possible way to qualitatively explain the behavior of the Phoenix CEL during the first 
half of its history is by mass-loss of not mixed enriched material as it is ejected from dying 
stars. As for the second half, a fast enrichment is possible if fresh, unprocessed gas is very 
scarce and stars are born from gas heavily contaminated by stellar ejecta. The step down 
in the SFR at ~ 6 Gyr ago could be related with the same gas shortage. T his qualitatiy e 



scena rio is compatible both with the fact that Phoenix has very little gas (lYoung et al 



20071 ). although it holds ~ 100 Gyr old stars. However, this is only a qualitative approach 



and more detailed models must be computed to better account for the process. 
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8. The SFH of Phoenix as a function of galactocentric radius 

The study of the spatial distribution of the stellar populations provides fundamental 
information about the formation and evolution of dwarf galaxies. Low surface brightness 
structures play a significant role in this issue. Being fossil records of the early stages of 
the evolution of the galaxies, they can provide key information about the composition and 



NOTE TO EDITOR: PLEASE, PUT FIGURE [lU] HERE 
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structure of dwarf galaxies when they started to form stars. In particular they can reveal if 
a real halo formed only by old, low metallicity stars is present (see 31]) ■ To this end, stellar 
age distribution must be accurately known in order to reach well-founded conclusions. This 
information is contained in the SFH ip{t,z). As we mentioned in ^ there are evidences 
that Phoenix has a stellar population gradient. In this section we will study such a gradient 
by solving ilj{t, function of the galactocentric radius. 

To derive the SFH as a function of radius, the observed stars were assigned to different 
regions of the gala xy. These regions were define d by fitting isophotes to ground based 



images of Phoenix (iMartmez-Delgado et al 



19991 ) obtained with the LCO 100-inch telescope 



(Chile). To fit the isophotes, images observed in I were combined and the resulting image 
was smoothed with a Gaussian filter of a = 50 pixels (~ 13"). Ellipses were fitted to the 
isophotes with ellipticity, position angle, and center as free parameters. A set of 6 ellipses 
of increasing semi-major axis were selected to divide the field in regions, each one including 
a statistically representative number of stars. 

Figure [TT] shows the selected ellipses overplotted on the ground based image of the 
galaxy. Figure [12] shows the six CMDs for the elliptical regions shown in Fig. [TTl Table 
[5] shows the fitted parameters of the selected ellipses and the number of stars within the 
corresponding elliptical region. For each ellipse (# 1 is the inner one), column 2 gives the 
semi-major axis in pc, column 3 gives the position angle, column 4 gives the ellipticity, 
and column 5 gives the number of stars between that ellipse and the former. Note the 
' 90° change in orienta t ion b etween ellipses #2 and #4, which was already noted by 



Martfnez-Delgado et al. 



(119991 ) and that suggests two different galactic components. 
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The SFHs for CMDs shown in Fig. [121 were solved as described in ^ In particular, the 
observational effects were simulated on the sCMD according to the completeness tests for 
each region. Solutions are shown in Fig. [131 Figure UM shows iplt) (left panel) and the CEL 
(right panel) for each elliptical region. The CEL plots stop at ages for which no significant 
star formation is found. Regions #1, #2, and #3 represent the main contribution to the 
global SFH. Regions #4 to #6 show mostly star formation at old and intermediate ages 
but the metallicities are similar to those of inner regions for the same age. 



E16 



E17 



Figure [15] illustrates the stellar population gradient in Phoenix. Upper, middle, and 
lower panels show the mean values of the SFR {< ip >), age, and metallicity (< [M/H] >) 
as a function of the galactocentric distance. Although a step down is observed in < ip > 
in the outermost region, its behaviour is compatible with a single exponential decrease, 
without any significant change of slope. < [M/H] > goes slightly down in the inner 200 pc 
and remains approximately constant from that distance on. 

Table [6] gives a summary of the main properties of the SFH of Phoenix as a function of 
galactocentric distance. The limits of the semi-major axis (in pc) for each elliptical region 
is given in column 1; columns 2 to 5 give the values of < ip >, < [M/H] >, and < age >, 
respectively. 
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Figure [16] shows the ages corresponding to the 90th-percentile (ego) and the 5th- 
percentile (es) of the age distribution, ego is the age such that 90% of the star formation has 
happened more recently than it, and similarly e^. In other words, ego and es are indicative 
of the age at which the star formation start and finishes, respectively. The flat distribution 
of ego shows that the oldest population has almost the same age all over the galaxy. The 
steep behaviour of shows that the younger the stars, the more concentrated they are. 
But, in relation to the existence or not of a real old halo in Phoenix, the most interesting 
characteristic of the behaviour is that it shows a smooth variation with no sharp step. 

In summary, all together, the smooth gradients or flatness oi < ip >, < [M/H] >, 
< age >, ego, and are difficult to reconcile with the presence of an halo in Phoenix, 
as a break in the stellar population gradients would be expected in such a case 



( iTikhonov. Galazutdinova. fc Drozdovsky 



20051 ). The presence of intermediate age stars in 



the external regions neither supports the presence of an old halo in the galaxy. The change 
in the orientation of the ellipses observed between regions #2 and #4 does not correlate 
with any sharp change in the slope of any of the mean and percentile values obtained 
here. The orientation of the inner ellipses could be then related to the distribution of the 
youngest stars in the last star formation event. Therefore, we can state that at a distance 
of about 500 pc from the center, no traces of an old halo are observed in Phoenix. 
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9. The evolution of the galaxy structure 



The behaviour of 65 (Fig. [16]) shows that the age of youngest stars decreases as we 
move to inner regions in the galaxy. This could be produced if stars form in an inner region 
and then migrate to outer ones. Alternatively, if stars are formed in situ, it is compatible 
with a shrinking star forming region centered in the galactic center and initially extending 
out to the galaxy outskirts. Contraction of this region would be the simple consequence 
of pressure support reduction as gas supply diminishes. As we have already mentioned, 



Stinson et al. 



(I2OO9I ) have shown on the basis of theoretical models that the latter would be 



the most efficient mechanism at work in dwarf galaxies. 

Our results for Phoenix allow us to study the way in which the galaxy evolves and, in 
the case that stars are formed in situ, to compute the rate at which the star forming region 
shrinks. The morphological evolution of the galaxy can be described in terms of what 
would be seen if it would be observed at different look-back times or red-shifts. To this 
purpose, the SFRs shown in Fig. [THhave been given as input to lAC-star to compute the 
corresponding synthetic populations as would be observed at different times. The surface 
mass density in alive stars as a function of radius and time, pm{T,t), has been derived using 
time steps of 1 Gyr. Results are shown in Figure [T7I for a sample of six values of look-back 
time. The most noticeable results from this figure are the reverse slope for the early epoch 
and the break at about 350 to 450 pc, also for the early epoch. A small depletion is also 



visible in the innermost ellipse at all 



[see e.g. 



Aparicio &: Tikhonov 



200G 



epochs, which 



Aparicio et al. 



las been fo und in other dwarf galaxies 



2000 



2OOII ). Exponential fits to data 



corresponding to ellipses 2 to 5 are also shown. Innermost and outermost ellipses have been 
excluded due to the mentioned departures from the intermediate ellipses trend. 

The scale length, am{t), and central density, Pm,o{t), both as a function of time, have 
been computed from fits similar to the ones shown in Figure [I7| but using all the available 
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time snapshots (from to 12 Gyr with 1 Gyr time step). Results are plotted in Figure 
TH] which, for more clarity, shows am{t) in two different scales as well as a~^(t), i.e. the 
slope of the lnpin(r, t)vs.r fit. A red-shift scale is also given according to the cosmological 



param eters for a flat Einstein-de Sitter universe and the 3-years WMAP data (jSpergel et al. 



2003). 



There are several interesting things to note from this figure. The most striking one, 
already seen in Figure [T7I is that the slope of the mass density distribution, a~^(t), was 
negative in the early epoch and changed sign at about 10.6 Gyr ago. In other words, mass 
density in stars was higher in outer regions than in inner ones during the early phase of 
the galaxy evolution. Besides that, the scale length of the alive stars component mass has 
steadily and smoothly decreased since 11 Gyr ago while the central density has increased 
since the earliest epoch. Quadratic fits can be done to a^^{t) and pm,o{r,t), but it is 
remarkable that linear fits with a break somewhere between 6 and 7 Gyr ago also provide 
a good representation of data. This epoch is coincident with the time at which chemical 
enrichment departs from a closed box (see Section [Tj). We have pointed to the possibility 
that fresh gas consumption would have been the cause of the latter. If shrinking is produced 
by the reduction of gas supply, it is not surprising that fresh gas consumption would 
produce some change in the shrinking rate. However, theoretical models are necessary to 
ascertain if this is the case. 

In summary, if stars form in situ in Phoenix, Figures [T7] and [18] indicate that the onset 
of star formation took place all over the galaxy, up to a distance of about 400 pc from the 
center at least, and with a relatively larger intensity in intermediate to outer regions, rather 
that in the innermost ones. After a few Gyr, star formation in inner regions started to 
dominate. Besides that, the star formation has become more and more concentrated, slowly 
shrinking down as time goes on. This scenario shows good agreement with model results by 
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Stinson et al.l (|2009| ) and seems to rule out the possibility that the galaxy holds a true old 
halo or, at least, that it could be significant at 500 pc from the center, equivalent to about 
2.5 scale- lengths. 



10. Summary and conclusions 

Photometry of resolved stars for the Phoenix dwarf galaxy using WFPC2@HST with 
the F555W and F814W filters have been presented. Color-magnitude diagrams, distance, 
mean metallicity, and the SFH, both global and as a function of radius has been obtained. 
This information allowed us to study the main features of the galaxy and the spatial 
distribution of its stellar populations. Results can be summarized as follows: 

(1) Using the TRGB, a distance modulus of (m — M)q = 23.09 ± 0.10, corresponding 
to d = 415 ± 19 kpc has been obtained. A mean metallicity [Fe/H] = —1.87 ± 0.06 dex, 
{Z = (2.6 ± 0.3) X 10^^) for the RGB stars has been also derived. 

(2) Using lAC-pop, the star formation rate as a function of time and metallicity, ilj{t,z), 
has been obtained, globally and as a function of radius. Phoenix has stars of all ages with 
a metallicity range of Z = 0.0001 — 0.002. Intermediate age and old stars can be found at 
all galactocentric distances. 

(3) Phoenix has had ongoing but gradually decreasing star formation over nearly a Hubble 
time. Star formation started in Phoenix at a very early epoch: 10% of the stars ever born 
in Phoenix were already formed about 13 Gyr ago. 50% of the total star formation had 
already taken place by 10.5 Gyr ago. Then, star formation continued at a significant level 
until 6 Gyr ago, forming an additional 35% of the stars, and at a very low level till around 
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100 Myr ago. 

(4) The distribution of mass density in alive stars and its evolution with time has been 
studied. If the stars we observe today were formed in situ in Phoenix, the onset of the star 
formation took place all over the galaxy, to a distance of at least 400 pc from the center. 
A relatively larger rate is observed in intermediate to outer regions in the earliest star 
formation episode, indicating that the star formation onset occurred preferentially in those 
external regions. 

(5) After a few Gyrs, the star formation in inner regions started to dominate, although it 
continued at a relatively high rate also in intermediate and outer regions for several Gyrs. 
At about 7 to 8 Gyrs ago, star formation stopped in the outermost regions (about 0.5 kpc 
from the center). From that point on, it was progressively stopping closer to the central 
regions as time went on, until the present, in which recent residual star formation activity is 
observed in the innermost region. As a result of this process, the scale-length of the stellar 
mass density distribution decreases with time. 

(6) The former sketch is compatible with a scenario in which the star formation region 
envelope slowly shrinks as time goes on, possibly as the natural result of pressure support 
reduction as gas supply diminishes. This results in a) old and intermediate- age stars present 
at all galactocentric distances but young stars only in the inner regions, b) a positive 
galactocentric gradient of age of the youngest population and a positive small gradient in 
the mean stellar age. All these results are in agreement with our observations. 

(7) The chemical enrichment law has been obtained from ijj{t,z). Chemical enrichment 
was very slow until about 6 Gyr ago, when metallicity increased abruptly. The rough 
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coincidence with the time at which the star forming region envelope start to shrink allows 
us to speculate that, from that time on, the fraction of fresh, unprocessed gas involved 
in subsequent star formation dramatically decreased. Detailed theoretical models are 
necessary to ascertain the plausibility of this suggestion. 

(8) No traces of a true, old halo are apparent in Phoenix either in its stellar age 
distribution or in the stellar mass density distribution, at least out to 0.5 kpc (about 2.5 
scale-lengths) from the center. The time evolution of the mass density distribution does not 
show traces of such structure either. 

We are very grateful to E. Skillman, P. B. Stetson, and K. McQuinn for their 
suggestions to improve this paper. The authors are funded by the lAC (grant P3/94) and 
by the Science and Technology Ministry of the Kingdom of Spain (grants AYA2004-3E4104 
and AYA2007-3E3507). 



-31 - 

REFERENCES 

Aparicio, A., Carrera, R, & Martinez- Delgado, D. 2001, A J, 122, 2524 

Aparicio, A., Dalcanton, J. J., Gallart, C, & Martinez- Delgado, D. 1997, AJ, 114, 1447 

Aparicio, A., & Gallart, C. 1995, AJ, 110, 2105 

Aparicio, A., & Gallart, C. 2004, AJ, 128, 1465 

Aparicio, A., & Hidalgo, S. L. 2009, A J, accepted. 

Aparicio, A., & Tikhonov, N. 2000, AJ, 119, 2183 

Aparicio, A., Tikhonov, N., & Karachentsev, I. 2000, AJ, 119, 177 

Baggett, S., et al. 2002, in HST WFPC2 Data Handbook, v. 4.0, ed. B. Mobasher, 
Baltimore, STScI 

Battinelli, P., & Demers, S. 2000, AJ, 120, 1801 

Battinelli, P., & Demers, S. 2006, A&A, 447, 473 

Battinelli, P., Demers, S., & Artigau, E. 2007, A&A, 466, 875 

Battinelli, P., Demers, S., & Kunkel, W. E. 2006, A&A, 451, 99 

Bellazzini, M., Ferraro, F. R., SoUima, A., Pancino, E., & Origlia, L. 2004, A&A, 424, 199 

Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F., & Nasi, E. 1994, A&AS, 106, 275 

Binggeli, B., Tarenghi, M., & Sandage, A. 1990, A&A, 228, 42 

Canterna, R., & Flower, P. J. 1977 ApJ, 212, L57 

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 



-32- 

Carignan, C, Demers, S., & Cote, S. 1991, ApJ, 381, L13 

Castelli, F., & Kurucz, R. L. 2004, New Grids of ATLAS9 Model Atmospheres lAU 
Symposium 210, A20 

Charbonneau, P. 1995, ApJS, 101, 309 

Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328 

Dekel, A., & Silk, J. 1986, ApJ, 303, 39 

Demers, S., & Battinelli, P. 2002, AJ, 123, 238 

Demers, S., Battinelli, P., & Artigau, E. 2006, A&A, 456, 905 

Dolphin, A. 1997, New Astronomy, 2, 397 

Dolphin, A. E. 2002, MNRAS, 332, 91 

Gallart, C, Aparicio, A., Freedman, W. L., Madore, B. F., Martinez- Delgado, D., & 
Stetson, P. B. 2004, AJ, 127, 1486 

Gallart, C, Stetson, P. B., Hardy, E., Pont, F., & Zinn, R. 2004, ApJ, 614, L109 

Gallart, C., Martinez- Delgado, D., Gomez-Flechoso, M. A., & Mateo, M. 2001, AJ, 121, 
2572 

Gallart, C., ZoccaU, M., & Aparicio, A. 2005, ARA&A, 43, 387 
Girardi, L., Bressan, A., BerteUi, G., & Chiosi, C. 2000, A&AS, 141, 371 
Harris, W. E. 1991, ARA&A, 29, 543 
Harris, J., & Zaritsky, D. 2001, ApJS, 136, 25 

Held, E. V. 2002, Scientific Drivers for ESQ Future VLT/VLTI Instrumentation, 178 



-33- 

Held, E. v., Saviane, I., & Momany, Y. 1999, A&A, 345, 747 
Hidalgo, S. L., Marin-Pranch, A., & Aparicio, A. 2003, AJ, 125, 1247 

Hill, R. J., Ferrarese, L., Stetson, P. B., Saha, A., Freedman, W. L., Graham, J. A., Hoessel, 
J. G., Han, M., Huchra, J., Hughes, S. M., lUingworth, G. D., Kelson, D., Kennicutt, 
R. C, Bresolin, F., Harding, P., Turner, A., Madore, B. F., Sakai, S., Silbermann, 
N. A., Mould, J. R. & Phelps, R., 1998, ApJ, 496, 648 

Hodge, P. 1989, ARA&A, 27, 139 

Holtzman, J. A., Smith, G. H., & Grillmair, C. 2000, AJ, 120, 3060 

Irwin, M., & Tolstoy, E. 2002, MNRAS, 336, 643 

Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 

Komiyama, Y, Okamura, S., Yagi, M., Furusawa, H., Doi, M. et al. 2003, ApJ, 590, L17 

Leisy, P., Corradi, R. L. M., Magrini, L., Greimel, R., Mampaso, A., & Dennefeld, M. 2005, 
A&A, 436, 437 

Letarte, B., Demers, S., BattineUi, P., & Kunkel, W. E. 2002, AJ, 123, 832 
Madore, B. F., & Freedman, W. L. 1995, AJ, 109, 1645 
Magrini, L., et al. 2003, A&A, 407, 51 

Martmez-Delgado, D., Gallart, C, & Aparicio, A. 1999, AJ, 118, 862 

Mayer, L., Governato, F., Colpi, M., Moore, B., Quinn, T., Wadsley, J., Stadel, J., & Lake, 
G. 2001, ApJ, 559, 754 

Mayer, L., Mastropietro, C., Wadsley, J., Stadel, J., & Moore, B. 2006, MNRAS, 369, 1021 



Menzies, J., Feast, M., Whitelock, P., Olivier, E., Matsunaga, N., & da Costa, G. 2008, 
MNRAS, 385, 1045 

Mighell, K. J. 1999, ApJ, 518, 380 

Minniti, D., & Zijlstra, A. A. 1996, ApJ, 467, L13 

Minniti, D., Zijlstra, A. A., & Alonso, M. V. 1999, AJ, 117, 881 

Moore, B., Calcaneo-Roldan, C, Stadel, J., Quinn, T., Lake, G., Ghigna, S., & Governato, 
F. 2001, Phys. Rev. D, 64, 063508 

Noel, N. E. D., & Gallart, C. 2007, ApJ, 665, L23 

Noel, N. E. D., Gallart, C., Costa, E., & Mendez, R. A. 2007, AJ, 133, 2037 
Ortolani, S., & Gratton, R. G. 1988, PASP, 100, 1405 

Peimbert, M., Colin, P., & Sarmiento, A. 1994, Violent Star Formation, from 30 Doradus 
to QSOs, 79 

Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 
Purcell, C. W., Bullock, J. S., Zentner, A. R. 2007, ApJ, 666, 20 
Sandage, A. 1971, ApJ, 166, 13 

Saviane, I., Rosenberg, A., Piotto, G., & Aparicio, A. 2000, A&A, 355, 966 
Schuster, H. E., & West, R. M. 1976, A&A, 49, 129 
Seth, A. C, Dalcanton, J. J., & de Jong, R. S. 2005a, A J, 129, 1331 
Seth, A. C, Dalcanton, J. J., & de Jong, R. S. 2005b, AJ, 130, 1574 



-35- 

Smith, G., Grillmair, C. J., & Holtzman, J. A. 1999, STScI, Proposal 6798 for HST 
Observations: "The Early Evolution of Local Group Dwarf Irregular Galaxies". 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Stetson, P. B. 1994, PASP, 106, 250 

Stinson, G. S., Dalcanton, J. J., Quinn, T., Gogarten, S. M., Kaufmann, T., Wadsley, J. 
2009, MNRAS, in press 

Tikhonov, N. A. 2005, Astronomy Reports, 49, 501 

Tikhonov, N. A. 2006, Astronomy Letters, 32, 149 

Tikhonov, N. A. 2006, Astronomy Reports, 50, 517 

Tikhonov, N. A., Galazutdinova, O. A., & Drozdovsky, I. O. 2005, A&A, 431, 127 
Turner, A., et al. 1998, ApJ, 505, 207 
Vansevicius, V., et al. 2004, ApJ, 611, L93 

van de Rydt, P., Demers, S., & Kunkel, W. E. 1991, AJ, 102, 130 
Young, L. M., & Lo, K. Y. 1997, ApJ, 490, 710 

Young, L. M., Skillman, E. D., Weisz, D. R., & Dolphin, A. E. 2007, ApJ, 659, 331 



This manuscript was prepared with the AAS macros v5.2. 



-36- 



3f: 

2- 

Or* 

-1- 
-2- 



-3- 



-4- 
-5-. 



J_L 




I I I 1 I V r I - I ^ 



-5 -4 -3 



-2 -1 
X 



1 



Fig. 1. — Image of Phoenix obtained with the 100-inch Du Pon^ telescope at Las Campanas 
Observatory (Chile). The Fields observed with HST are overplotted on the image. The 
center of the figure is located at the optical center of the galaxy. At the distance of Phoenix, 
1' corresponds to ~ 122 pc. North is up, East is to the left. 
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Fig. 2. — Uncertainties provided by ALLFRAME as a function of magnitude for F555W 
(top panel) and F814W (bottom panel) filters. 
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Fig. 3. — CMD of Phoenix. Absolute magnitudes are shown in the right axis. Main features 
of the CMD are shown on the right panel. 
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Fig. 4. — The luminosity function of stars redder than {V — I) = 0.8 (dashed hne) and its 
convolution with a Sobel filter of kernel [1,2,0,-2,-1] (solid black line). A vertical dash dotted 
line shows the convolution peak, placed at / = 19.14 ± 0.02 mag. Arrows indicate the peak 
half height. The error of the peak position is estimated as the width at half of height of the 
peak. 
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Fig. 5. — Globular cluster can didates "a", "b", and "c" (from top to bottom) proposed by 



Martfnez-Delgado et al. 



(119991 ). All of them seem to be background galaxies. 
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Fig. 6. — Observed CMD (red) overplotted on the CMD of the recovered artificial star CMD 
(black). Completeness is given on the right Y-axis 
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Fig. 7. — Regions of the CMDs used to obtain the SFH. The observed CMD is shown on the 
left and the sCMD, on the right. 
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Fig. 8. — The SFH ilj{t, z) of Phoenix. It contains a binary fraction of 40%. The height of the 
curved surface over the age-metalhcity plane gives z). The SFH depending only on time 
(ipit)) and metallicity {iIj{z)), are shown in the right and left vertical pannels, respectively. 
The CEL is the projection onto the horizontal plane. 
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Fig. 9. — Upper panel shows the SFH of Phoenix as a function of time only, ip{t). It is 
the same shown in the right-hand vertical panel of Fig. [HI Thin lines represent the rms 
dispersion of the solution. Middle panel shows the CEL in a logarithmic scale, obtained 
averaging iplt, z) in z for each value of age. In this case, thin lines represent the interval 
containing 68% of the distribution. Vertical dot-dashed lines show the times by which 10%, 
50% and 95% of the mass has been converted into stars. Bottom panel shows the normalized 
cumulative mass converted into stars. 
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Fig. 10. — The CEL of Phoenix (thick, black hne; thin black lines show the dispersion 
interval, see Figure [9]) compared with the results of several metallicity law models. Closed- 
box (dashed lines) for several final gas fractions /i/ are shown. Outflow models (dotted 
lines) with A = 1 are shown for the same fif used for the former. Infall models (dash- 
dotted lines) with a = 1 are also shown for fif = 0.1 and fif = 1. Only close-box and 
outflow models with Uf = 10~^ 10~^ are compatible with the CEL of Phoenix. (See 



Peimbert. Colin, fc Sarmiento 



19941 . for the definition of a and A parameters) 
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Fig. 11. — Ellipses selected to divide the galaxy field in regions for the study of stellar 
populations gradients on the image. The WFPC2 fields are also shown 
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Fig. 12. — CMD of Phoenix for each elhpse shown in Fig. [TTl Isochrones of t = 14 Gyr, z = 
0.0001 (red); t = 8 Gyr, z = 0.004 (green); t = 5 Gyr, z = 0.004 (blue); and t = 0.2 Gyr, z = 
0.0001 (cyan) have been overplotted. Labels 1 to 6 refer to the elliptical regions shown in 
Fig. [m radius increasing from 1 to 6. 
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Fig. 13. — iplt, z) for the elliptical regions shown if Fig. [TTJ The radius increases from panel 
#1 to #6 (See Table E]). 
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look-back time (Gyr) 

Fig. 14. — SFR (left panel) and CEL (right panel) for the six elliptical regions defined in 
Figure [Til Regions are labeled as in Figures [12] and [131 A boxcar average of width 1 Gyr 
and step of 0.1 Gyr has been used to smooth ilj{t) and z{t). In the CEL panel, metallicity 
law ends where no significant SFR is found. 
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Fig. 15. — Upper, middle and bottom panels show respectively the mean ip{t), mean age and 
mean metallicity of the stars as a function of galactocentric radius. In all the cases, vertical 
bars are rms dispersions and horizontal bars show in the interval of age embraced by the 
associated value. 
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Fig. 16. — The 90-percentile (ego) (open squares) and 5-percentile (65) (open circles) of the 
star age distribution. In both cases, vertical bars are rms dispersions and horizontal bars 
show the interval of age embraced by the associated value. 
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Fig. 17. — Distributions of surface mass density in alive stars as a function of galactocentric 
radius for several look-back times in Phoenix. Lines show fits to data for ellipses 2 to 5. 
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Fig. 18. — Central surface density and scale length of the distribution of surface mass density 
in alive stars as a function of look-back time in Phoenix. Panel (A) shows the central density; 
open circles represent the values obtained from the fit to the alive star mass distribution, 
while asterisks provide the average values directly measured in the inner ellipse (see Fig. 
[TT]) . Red line shows two linear fits performed to the values corresponding to look-back times 
larger and smaller than 6.5 Gyr, respectively. Blue lines represent a quadratic fit. Panel 
(B) shows the scale length. Panel (C) shows the same with a different vertical scale, for 
more detail. For illustrative purposes. Panel (D) shows the inverse of the scale length, to 
emphasize the fact that the slope of the alive star mass distribution changes sign at about 
10.6 Gyr ago. In panels (B) to (D), red lines show two linear fits to performed to the 
values corresponding to look-back times larger and smaller than 6.5 Gyr, respectively, while 
blue lines represent a quadratic fit, also to a^. 
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Table 1. Main properties of Phoenix dwarf galaxy. 



Property Value Refs. 



Galactic coordinates 6) 272? 19, -68? 95 5 

Equatorial coordinates (02000, feooo) 0l''51™03!3, -44°27'11" 5 

Distance modulus (m - M)o 23.1 ± 0.1, 23.0 ± 0.1, 23.21 ± 0.08, 23.11, 23.09 ± 0.10 6,4,2,3,7 

Reddening AyAl 0.062, 0.036 1 

Absolute Magnitude M°, -9.8 ± 0.4 4 

Total mass (10® M©) 4.4 ± 0.3 7 

[Fe/H] -1.81 ±0.10 2,3 

Average star formation rate {10~^ MQyr~^pc~^) 0.99 ± 0.08 7 



References. — (1) Cardelli, Clayton, & Mathis 1989; (2) Held, Saviane, & Momany 1999; (3) Holtzman, Smith, & 
Grillmair 2000; (4) Martmez-Delgado, Gallart, & Aparicio 1999; (5) Schuster & West 1976; (6) van de Rydt, Demers, & 
Kunkel 1991; (7) This work. 
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Table 2. Observation Log. 



Visit 


Field 


Cycle 


Proposal ID 


Number of orbits 


Filter 


Number of images 


Exp. time (s) 


2 


inner 


6 


6798 


2 


F555W 


6 


4900 


2 


inner 


6 


6798 


2 


F814W 


6 


4900 


1 


inner 


9 


8796 


3 


F555W 


6 


7200 


2 


inner 


9 


8796 


4 


F814W 


8 


9600 


3 


outer 


9 


8796 


5 


F555W 


11 


11900 


4 


outer 


9 


8796 


6 


F814W 


13 


14300 
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Table 3. Photometric calibration terms. 



# CCD 


Z 


C 


K X 10-* 




PC (V) 


21.623 





1.27 


0.012 


PC (I) 


20.829 


-0.053 


0.94 


0.004 


WFC2 (V) 


21.740 


-0.032 





0.003 


WFC2 (I) 


20.888 


-0.023 





0.004 


WFC3 (V) 


21.796 


-0.063 





0.003 


WFC3 (I) 


20.896 


-0.033 





0.004 


WFC4 (V) 


21.657 





2.19 


0.003 


WFC4 (I) 


20.840 


-0.030 


1.37 


0.004 
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Table 4. Mean and integrated values derived from il^it^z). 



SEL 


Total Mass 


< V > 


< [M/H] > 


< age > 




(106 Mq) 


(lO-^Moj/r-- V"^) 




(Gyr) 


BaSTI 


4.4 ±0.3 


0.99 ± 0.08 


-1.7±0.2 


9.4 ± 1.1 


Padua 


4.8 ±0.4 


1.04 ±0.07 


-1.7±0.1 


10.0 ± 1.2 
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Table 5. Fitted parameters for the ellipses in Fig. [TTl 



ellipse # 


sma (pc) 




e 




1 


104.2 


-92 


0.5 


8640 


2 


151.8 


-80 


0.2 


12235 


3 


199.0 


-24 


0.1 


5578 


4 


261.4 


-9 


0.3 


3011 


5 


364.8 


-9 


0.3 


3013 


6 


573.2 


-9 


0.3 


1585 
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Table 6. Mean values as a function of radius. 



sma 


<ip> 


< [M/H] > 


< age > 


(pc) 


(lO-^Mgyr- V"^) 




(Gyr) 


- 104.2 


1.61 ±0.09 


-1.57 ±0.13 


8.2 ±1.3 


104.2 - 151.8 


1.42 ±0.09 


-1.62 ±0.14 


8.8 ± 1.2 


151.8 - 188.0 


1.26 ±0.09 


-1.66 ±0.13 


9.6 ± 1.1 


188.0 - 261.4 


1.08 ±0.17 


-1.74 ±0.11 


10.7 ±0.9 


261.4 - 364.8 


0.64 ± 0.06 


-1.74 ±0.13 


11.2 ±0.7 


364.8 - 573.2 


0.30 ± 0.04 


-1.80 ±0.08 


10.8 ±0.7 



